Thermally stimulated H emission and diffusion in hydrogenated amorphous silicon 
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We report first principles ab initio density functional calculations of hydrogen dynamics in hy- 
drogenated amorphous silicon. Thermal motion of the host Si atoms drives H diffusion, as we 
demonstrate by direct simulation and explain with simple models. Si-Si bond centers and Si ring 
centers are local energy minima as expected. We also describe a new mechanism for breaking Si-H 
bonds to release free atomic H into the network: a fluctuation bond center detachment (FBCD) 
assisted diffusion. H dynamics in a-Si:H is dominated by structural fluctuations intrinsic to the 
amorphous phase not present in the crystal. 

PACS numbers: 66.30. Jt,66.30.-h,61.43.Dq 



I. INTRODUCTION 

Hydrogenated amorphous silicon is one of the most im- 
portant electronic materials, and is used in applications 
ranging from TFTs in laptop displays to solar photo- 
voltaics and IR imaging/detection. Hydrogen dynam- 
ics is key to creation and annihilation of defects, and is 
also linked to light induced degradationii^i^i^. Although 
there have been several studies, a complete picture of H 
dynamics is not yet available. 

A widely held picture posits that unbonded H hops 
among various attractive sites before capture at a dan- 
gling bond^i^iii^i^iiS. Au intermediate low-energy path- 
way involving a metastable dihydride structure has also 
been reportedii. However, by calculating hopping rates 
for different trap sites, Fedders argued that thermal mo- 
tion of hydrogen does not proceed from dangling bond to 
dangling bond via bond center (BC) sites and showed dif- 
fusion through intermediate levels to be insignificantji^. 
In addition, Su et afl^ proposed that Si-H bonds do not 
spontaneously release H, but rather require the mediation 
of an external agency: in their case a five-fold or "float- 
ing" bond. In crystalline Si, the importance of lattice dy- 
namic activated diffusion has been reportedi^ii^. Buda 
et ali^ has shown diffusion of H in the form of jumps 
from BC site to another BC via intermediate hexagonal 
or tetrahedral sites. 

Previous work on a-Si and a-Si:H showed that the net- 
work dynamics is in some ways quite different from the 
crystal. In particular, the disorder of the network allows 
fluctuations in the positions of atoms leading to the inter- 
esting observation of "coordination fluctuation" . It has 
been observed in an early first principles simulation that 
even at T=300K, the number of floating (fivefold) bonds 
fluctuated between zero and ten in a 216 atom cell, in a 
1.8 ps simulations^. This work has been updated, and 
similar effects have been observed in networks including 
H. It was also found that most of the atoms in the lattive 
eventually participated in these fluctuationsi^. 

In this Letter, we report an ab initio simulation which 
reveals the key role of thermal motion of Si atoms in 
driving H diffusion. We have undertaken accurate simu- 
lations including static lattice simulation (in which Si 



atoms were frozen) and extended thermal simulation. 
The static Si lattice simulations shows no H diffusion 
as compared with the dynamic lattice case, suggesting 
that the motion of the "Si-sublattice" is important to 
the H dynamics. A key feature of our work is that we de- 
termine diffusion mechanisms directly from thermal MD 
simulation, not by imposing a conventional hopping pic- 
ture among wells (traps) with varying depths. The prin- 
cipal conclusion of the paper is that the dynamic lattice 
(particularly the motion of pairs or triplets of Si atoms 
with a BC H present) is a primary means for ejecting 
atomic H into the network. This mechanism could not 
be easily inferred from phenomenological kinetic equation 
models of H transport^, though it should readily emerge 
using a method devised to discover rare (long time scale) 
events like the Activation-Relaxation Technique imple- 
mented with ab initio interactionsiS. 

We determine an essential mechanism for H diffusion 
in the dynamic lattice, namely, "Fluctuating Bond Cen- 
ter Detachment" (FBCD) assisted diffusion: if the H is 
initially covalently bonded to a Si atom, it stays bonded 
with it until another Si comes in the vicinity and makes 
an instantaneous or fluctuating BC configuration. This 
event is followed by a switching of H from the covalent 
bond to the new Si to either form another Si-H bond or 
hop, depending upon the local environment. This process 
is important both as a means for the network to gener- 
ate free H and to create dangling bonds. The mechanism 
we report here is undoubtedly not the sole means of ob- 
taining H diffusion, but is predominant in accurate and 
relatively extended MD simulations. 



II. METHODS 

To understand the diffusion of H at different tempera- 
tures we have used a 71 atom model aSieiHiQ. We gen- 
erated this model by removing three Si atoms from a 
well relaxed defect free aSiQ4 mode^. We then termi- 
nated all the dangling bonds except two with hydrogen. 
The newly formed supercell is then relaxed using conju- 
gate gradient to form a aSieiHio model (61 Si and 10 H 
atoms). The computed Si-H partial correlation has the 
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FIG. 1: Si-Si and Si-H partial pair correlation for aSiaiHio 
and Si-Si partial pair correlation for 05164 models. 



In our simulations, H in a-Si:H is either passivating 
Si dangling bonds, or occupying other sites, such as BC 
sites (the most common) or certain other locations that 
provide energy minima for H. To provide some informa- 
tion about the dynamics of our cells, we have considered 
dynamics near BC conformations in aSieiHio. In Fig- 
ured we have plotted, p{R, 9) oc Y.^ Ej S{R-Ri{t))5{e- 
6j{t)), a distribution function indicating time spent in 
different parts of the R — 6 configuration space, where 
Ri{t) is the distance between two Si atoms at a time, 
and 9j{t) is the (Si-H-Si) angle formed by two silicons 
and a central H at a given time. This result shows pre- 
ferred values of R, ranging from 3.0 - 4.0 A and 9 (Si- 
H-Si bond angle) in the range of 110° - 180° where the 
H atoms spends most time. The hydrogen trapping time 
is highest on the red and lowest on the green regions. 
Evidently, the energy surface is rather weakly dependent 
upon 9 over a wide range of angles. 



expected peak at 1.5 A and the Si-Si partial pair cor- 
relation show a first peak at 2.35 A with a subsequent 
minimum which remains unchanged with hydrogenation 
as shown in Fig. [T] The newly generated model retains 
a tetrahedral structure and fourfold coordination (for Si) 
as in the original. We have also considered additional 
models, a 138 atom model aSii2oHi^ and a 223 atom 
model aSi2i4H^^^, and studied H dynamics there to as- 
sure ourselves that finite-size effects were unimportant. 

The simulations were performed using SIESTA^i in 
the generalized gradient approximation for the ex- 
change (GGA) using a parametrization of Perdew, Berke, 
and Ernzerho:^. Norm conserving TrouUier-Martins 
pseudopotentiala^^ factorized in the Kleinman-Bylander 
form^^ were used. We used a single C polarized basis set 
for Si valence electrons and double C polarized basis for H 
allowing sufficient flexibility in the basis set. We solved 
the self-consistent Kohn-Sham equations by direct diag- 
onalization of the Hamiltonian and a conventional mix- 
ing scheme. Every structure in this report was relaxed 
using conjugate gradient (CG) coordinate optimization 
until the forces on each atom were less than 0.02 eV/A. 
The r point was used to sample the Brillouin zone in all 
calculations. 



III. SIMULATIONS 

Our strategy has been to study the dynamics of H 
atoms in small cells which appear to represent the local 
topology of the amorphous network, using highly accu- 
rate methods. The use of such methods precludes simu- 
lations much exceeding several picoseconds, particularly 
because of the small mass of H and the commensurately 
small time step needed to integrate the equations of mo- 
tion. 




2 2.5 3 3.5 4 4.5 5 

Si-Si distance (A) 



FIG. 2: Dynamics of H near BC conformations. The contours 
give normalized temporal distribution of Rsi-Si as a function 
of Si-H-Si bond angle from aSieiHio for a total of 6.25 ps 
at T=1000K: red, most frequent visitation and green least 
frequent. The plot is averaged over all H atoms in the cell 
when in a BC conformation. 



To analyze the role of thermally induced Si motion in 
driving H diffusion, we considered 6.25 ps simulations (25 
000 MD steps each with time step r=0.25 fs) for aSieiHio 
at T=300K and T=1000K. The difference in the H diffu- 
sion for these two cases is easily extracted from the time- 
averaged mean squared displacement of H, which is de- 
fined as (a2)t,„e = E.=1 m) - rl(0)p, 

where Nm d is the total number of MD step, Nh and (i) 
are total number and coordinates of H at time t respec- 
tively. This yielded an average mean square displacement 
of 0.14 A^ for T=300K and 1.06 A^ for T=1000K. Simi- 
lar analysis on Si gives an average mean square displace- 
ment of 0.06 A2 and 0.57 A^ for T=300K and T=1000K 
respectively. 
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IV. DISCUSSION 

In the case of aSigiHio at T=1000K, we have ob- 
served 9 bond breaking events which are accompanied 
by rapid bond switching in the 6.25 ps simulation time. 
For aSi2iiHg at BOOK we have observed 5 bond break- 
ing events in the 10 ps simulation time; all are FBCD 
assisted. Only one bond breaking event we observed for 
aSiQiHiQ at BOOK for the total simulation time of 6.25 ps. 
We observed no bond breaking event for aSii2oHig model 
at 300K, 500K and TOOK temperatures. After the bond 
breaking, all of the events lead to the hydrogen passivat- 
ing a dangling bond. To discuss the FBCD mechanism 
in detail, we have selected two H from aSieiHiQ model 
namely Hgy and Hgg at T=1000K, which diffuse through 
via FBCD. To analyze the role of the thermal motion of 
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FIG. 3: Her switching the bond from Siss to Sisi assisted 
by the close approach of Sias. a) shows the swhiching of H 
bond between the two Si. b) and c) shows the fluctuation 
of distances of Siaa and Sisi from Sias and variation in the 
corresponding angles centered at He? respectively. 

In Figure [HI we show a situation in which Hgr initially 
bonded to Sias switched to Sisi. This event follows the 
close approach of Sisi and Sias to the Sias-Hgr bond and 
forms a fluctuating bond center conformation. As the 
distance Siss-Sisi and Sias-Siaa change from 4.8 A to 
3.1 A and from 5.0 A to 4.0 A respectively as shown in 
Fig-EIb), the bond angle Siss-Hey-Sisi changes from 140° 



to 80°. As the same time, the bond angle Siss-Hgr-Sisa 
also changes from 160° to 95°. These introduce change 
from a BC configuration, which is a local energy minima, 
to energetically unstable part of the R — configuration 
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FIG. 4: Another H Bond switching event, a) shows Hes 
switching bond from Si48 to Siae. b) and c) shows the fluctu- 
ation in the distances of Sise and Siag from Si48 and variation 
in the corresponding Si-H-Si angles respectively. 

A second case for the FBCD assisted diffusion in- 
volves Hgg similar to the previous example. As shown in 
Fig. Hl^a) , we observed a bond breaking process in Hgs , 
a situation where Hgs which was initially bonded to Si48 
and later switched to Siag. The switching occurs due 
to the close approach of Siag to the Si48-Hg8 bond and 
forms a fluctuating bond center conformation. As the 
distance Si48-Sia9 changes from 4.2 A to 3.0 A as shown 
in Fig. m^b), the bond angle Si48-Hg8-Sia9 changes from 
a bond center configuration 150° to 89°, compelling the 
H to diffuse and form a bond with another Si. 

In Fig. [51 we show the trajectory of Hg8 in FBCD 
assisted diffusion. The average positions of IIg8 are shown 
in yellow, while Si48, Siag, and Siag are shown with blue, 
orange, and green colors respectively at three different 
time steps of the simulation. IIg8 is initially bonded to 
Si48, becomes a fluctuating bond center H with approach 
of Siag, is ejected and eventually bonds to Siag. 

The FBCD mechanism is reminiscent of that of Su et 
aZ.^^, also depending upon the intercession of a Si not 
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in which no Si was overcoordinated. Ejection of H is a 
more subtle process than just changing coordination, and 
depends on the the local geometry (R, 9). Finally, all 
such FBCD conformations arise from fluctuations, and 
are thus short lived^®. The mechanism of Su et al is a 
special case of FBCD. 

In conclusion, we have demonstrated the nature of H 
diffusion in a-Si:H by direct simulation and with the aid 
of a model to develop a fairly simple picture of H motion. 
H emission is stimulated by Si motion, and the FBCD 
mechanism is shown to be important both for stripping 
off H chemically bonded to Si (thus creating "free" atomic 
H), and of course for creating Si dangling bonds. Our 
work is consistent with analogous studies in c-Si^^, and 
is a generalization of the work of Su and Pantiledes^^. 
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FIG. 5: Fluctuating Bond Center Detachment and diffusion. 
Trajectory (averaged over 400 steps) of Has are shown in yel- 
low, while Si48, Siae, and Sisg are shown with blue, orange, 
and green colors respectively for three different time steps of 
the simulation. 



part of the initial conformation. However, the FBCD is 
a more general process that may or may not increase the 
Si coordination. We found many cases of H detachment 
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